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ABSTRACT 


' The purpose of this work was to i 

develop a suitable computer program to 
correct computer compatible tapes, 
obtained from LANDSAT HSS data, and to 
study 3ts effect on percentage of correct 
classification The LOTORAK-3 program, 
developed by the Air Force Geophysics 
Laboratories and based on a semi-empirical 
model, was used ro calculate the 
atmospheric transmittance, in conjunction 
v'lth Turner's mode3 for correction of 
satellite data for atmospheric 
interference It improved the contrast 
between different natural targets in the 
MSS LandSAT data of Brasilia, Brazil It 
also improved the classif ication accuracy 
of sugar canes by about 9%, in the MSS 
LAIIDSAT data of Ribeirao Preto, Sao Paulo, 
Brazil. 


I I -!'"RODhCTION 

Radiation from remotely sensed 
objects in the earth's environment is 
attenuated in its pe.,3age through the 
atmosphere tilth the advent of LAKDSAT, 
Skylab and other ailvanced eartn monitoring 
satellites, it has fcacoma increasingly 
important to know to what ertent the 
atmosphere affects the classification 
accuracy, using satellite multispectral 
scanner (riSS) data, of earth resources 

The Brazilian Instj tute of Space 
Research (INPE) operates a data receiving 
and a data processing station for the 
LANDSAT satellite and processes flSS data 
in the form of photographs, as well as 
computer compatible tapes. The objectives 
of the present project v?ere* (1) To 
develop a suitable, computationally 
efficient, computer program to correct the 
tapes of uAWDSAT's multispectral scanner 
(MSS) data for atmospheric interference, 
and (2) , to study the effect of 
atmospheric interforonce on (A) the 


REPKODUCmiLITY OF THE 
OBIGINAL PAGE IS POOR 

quality of LAKDSAT MSS images for visual 
phoro-interprotatj on and (B) the 
percentage of correct classif icars on of 
the same data, using automatic pattern 
recognition techniques 


IT DETEajlINATIO\ OF 
ATMOSPHERIC TRANSMITTANCE 

The calculation of atmospheric 
transmittance is normally difficult, 
because of the complicated nature of the 
radiative transfer equatron and of the, 
usually incompletely knov/n, constituency 
of the atmosphere at tne time for winch 
calculations are to be made There have 
been innumerable investigations from v.hicn 
methods are derived for cajculeting 
atmospheric transmittance^ (Ij Line-by- 
line calculations- the foundation for the 
' line-by-line calculation is a ser of 
parameters that describe the molecular 
lines in vjhich radiation 3 s ab.,orbed and 
emitted Many techniques', developed in 
tliese calculations, are designed to 
minimize the computational effort Fven 
with the time-saving approximations 
applied to rnc line-hy-line cajculatiop, 
the costs are often excessive ror an 
accuracy beyond that required for 
correcting landsat mss data for 
atmospheric interference (2) Band-Model 
Calculations Although, in general, these 
techniques aie more efficient 
computationally than line-by-line 
caJculations, they still take considei ably 
more computer time than , for e <amplc , the 
LOKTRAN 3 method described below. (3) The 
Tiggregat e mctlod, developed by the 
Environmental Research Institute of 
Michigan (EFIM) , is .simply a compilation 
of tne special forms of the band-models, 
assembled in a program In this method, m 
certajn wavelength regions, the 
conventional hiodels are not used, but 
empirical fits to experimental data are 
used to detei mine the functional form of 
transmittance' Since its spectral range 
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Cl to 30 ym) does not include the visible 
wavelength region, it is not suitable for 
oir purpose. (4) LOWTRAN 3 Method The 
LOCTRAN 3 program^ is based on a semi- 
empirical model and calculates rhe 
transmittance (uveraged over a 20 cm 
interval) for a given atmospheric path^, 
at steps of 5 cm”'', from 350 cm~’ to 
40,000 cra-MO.25 to 28.5 vm) 

‘ The main assumptions made in this 
program arc that the atmosphere can be 
represented by a 33-layer model, and that 
the average transmittance over a 20 om~‘ 
interval (due to molecular absorption) can 
be represented by a sing3e parameter 
model, deuermined empirically using both 
laboratory transmittance data and. 
available molecular line constants. The 
absorption coefficients for water vapor, 
ozone, and the combined effects of the 
uniformly mixed gases (CO 2 , S 7 O, CHm , CO, 
1 ) 2 , and O 2 ) are digitized in this program 

A choice of si.^ model atmospheres 
[1952 U.S. Standard Atmosphere, Tropical 
(15°K) , Midlatitude Summer [45°N, Cfuly) , 
Midlatitude Winter (45°, January), 

Subarctic Sunmer (60°N, July), and 
Sub'arctic Winter (G0°K, January) , is 
given, with an option for a seventh model, 
which can be inserted as a set of 
radiosonde data. The aerosol attenuation 
IS calculated for a given visual range, 
using an interpolation [ extrapolation 
scheme, with two aerosol models, based on 
measurements of continental aerosols under 
moderate visibility conditions (5 km and 
23 km au sea level). The altitude, 
pressure, temperatire, vrater vapor density, 
and ozone density loi six model 
a-emospheres , as we] 1 a.s the number of 
particles per cin^ for tv;o haze mode]s 
(visual ranges of S and 23 )cra} , are 
provided as basic jnput data 

This program used for^ calculating 
the transmittance of the atmosphere 
because it is reasonably accurate, user- 
oriented, computationally very efficient, 
well documented and revised by the 
Air porce Geophy.sics Laboratories at 
reguJar periods of time, utiDizing recent 
laboratory measuieraents and ■cheorcstic.al 
calcu] ations. 


in CORRECT COM OF LliHDSAT MSS DATA FOR 
ATMOSPHERIC INTDRrERENCE 

Computer compatible tapes of LANDSAT 
,ISS data ijerc corrected using the 
'atmospheric modal for correction of 
spacecraft data' , as given by Turner 
(1972)'* This model assumes a plane- 
parallel atmosphere but can be applied to 
the case of a icalistic spherical 


atmosphere, if the nadir view ang) c is 
small, or if the satellite's altitude is 
not too great The target was assumed to 
be diffuse. The eight parameters needed by 
the model are v/avelength, azimuth and 
zenith angles of tne sun, azimuth and 
nadir view angles of the sensing device, 
he: gnt of the sensing device above the 
terrain, visual range at the ground, and 
average background terrain albedo . 

; The day, month, year, and local time 
of the satellite, along t’lth the longitude 
and latitude of the ground scene, were 
used in a separate computer program to 
obtain the solar azimuth and zenith 
angles The background terrain spectral 
albedo \/as estimated by the values of 
spectral reflectance of the known types of 
ground covers available in the literature 
It v;as found that changing the scanning 
angle of the satellite from 0 ° to 10 ° had 
almost no effect on the calculated 
transmittance of the atmosphere. Tl,e 
elapsed time on a B-6700 computer to 
correct 4 images, each of 512 > 512 pixels 
U e 1,048,576 pi/elsi for atmospheric 
^interference, vras about IS ininuces 
, ^ 

IV. EXPEPII4EHTAL RESULTS ^ 

LANDSAT MSS data of Biasilia, Brazil, 
vrere corrected Using this algorithm and 
available radiosonde data Figures I to 3 
'slio^z the original image in band 4 (0 5 to 
0 6 pm) , band 7 (0 8 to 1 1 ym] and a 
color composite, respectively. Figures 4 
to G shovi the respective figures 1 to 3 
correcced for atmospheric interference 

Mote that the corrected iiiuges have 
better contrast betvieen natural targets, 
thus making it easier to discriminate 
bcpween them by photo-inteipretation. In 
addition, one can see textural featur'=s of 
certain areas like vageuation, roads, etc 
better in the corrected images. 

To arrive at an estimate of chance 1.1 
percentage of 'correct classif ication, 
after correction for atmospheric 
interference, LAMDSAT MSS Data of RiLeirao 
Preto (July 1, 1977), was corrected for 
atmospheric interf orciico using this 
algorithm 'The DeparUnenl of Remote 
Sensing of Earth Resources of INPC had 
classified the sugar cancs in this area, 
using aerial photography The o: iginal end 
corrected images gave classification 
accuracies of sugar cancs to be 66 IS and 
77 3% respectively', using the single-cell 
option of image-.l .00 ns coniiiarcd to the 
results obtained by aerial photography 
Image-100 is a data processing system 
marketed by the General Electric Co to 
ext:-act thematic information and enhance 
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multispectral iriagery. Further erperrmenls 
need to be done to establish if correcting 
LAMDSAT HSS data for atmospheric 
interference significantly improves, 
statistically, the percentage of correct 
classification. Future plans include more 
experiments with the current algorithm, 
using a pixel-by-pixel maximum likelihood 
gaussian classifier, as well as a sample 
classifier, making the current algorithm 
computationally more efficient, and 
comparing the results to those of other 
algorithms of correcting satellite data 
for atmospheric interfeience reported in 
the literature, with respect to their 
effect on classification accuracy. 

We gratefully acknowledge the 
assistance of Dr Celso de Renna e Souza, 
Head, Informatics Division (INPE} and tne 
Department of Remote Sensing (INPE) for 
their assistance with this work. In 
addition, the authors would like to thank 
Dr. Nelson de Jesus Parada, Director of 
INPE, for his permission to publish ths s 
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Figure 1 . Original Image in Band 4 . 


Figure 2 - Original Image in Band 7 . 
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figure 3. Original Image - Composite color Figure 4. Corrected Image ift Band 4, 
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Pigure 5. Corrected Image in Baad 7. 
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Figure 6. Corrected Image Composite Color. 
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